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D . Abstract 

We analyse the azimuthal decorrelation of 'Mueller-Navelet' dijets produced in the 
pp collisions at Tevatron energies using a BFKL framework which incorporates dominant 
subleading effects. We show that these effects significantly reduce the decorrelation yet 
they are still insufficient to give satisfactory description of experimental data. However a 
good description of the data is obtained after incorporating within formalism the effective 
rapidity defined by Del Duca and Schmidt. 



An ideal test of log(l/x) dynamics is the "Mueller-Navelet" process [0,0,0] in hadron-hadron 
collisions with one backward jet separated from a forward jet by a large rapidity interval. The 
test has been under study for several years, but has recently been sharpened by the consid- 
erable improvement in the precision of the data Q. The relevant hard scale is specified by 
the transverse momenta of the jets (kijkz). An advantage of this process is that these mo- 
menta can be of comparable magnitude so the intervening DGLAP evolution is suppressed. 
The behaviour of the cross-section is therefore sensitive to the diffusion of transverse momenta 
of the intermediate gluons, which is a property of leading log(l/a;) dynamics [§, [|. The char- 
acteristic features of this process are (i) a rapid increase of the 2-jet production cross section 
with increasing incident energy and (ii) a weakening of the azimuthal back-to-back correlation 
between the jets as the rapidity interval is increased. The analyses have so far been within the 
leading order log(l/x) BFKL framework, and indicate much more decorrelation than shown 
by the data. The improvement of the 2-jet data has been accompanied by progress in under- 
standing and quantifying the effect of sub-leading log(l/x) corrections. The full NLO BFKL 
(log(l/x)) contributions have been obtained and found to be large 0, ||. However the ma- 
jor contribution can be resummed to all orders in a straightforward way |IUL O, O. To be 
precise the resummation is effected by requiring that the virtuality of the intermediate gluons 
is dominated by their transverse momenta squared, which although formally sub-leading, is 
implicit in the derivation of the BFKL equation. We call this the consistency condition (CC). 
The resummation is found to stabilise the solution of the NLO BFKL equation. The purpose 
of this paper is to incorporate those subleading effects generated by the consistency constraint 
in the theoretical description of the dijet production in hadronic collisions. 

The cross-section describing the production of a forward jet (xi,kf) and a backward jet 
(x2, k%), with relative azimuthal angle — vr, may be written 



da 



da 



dx\d%2dk\dk\d<\) X\X 2 \ Xl ^ cS ^ Xl ' i) X2 / eff ( a; 2? ^2) ^2^2^' ( ) 

where / c fr = + §(<? + <?) is the effective parton distribution of the incoming proton (antiproton). 
Note, that with this convention the back-to-back configuration of the jets in the transverse 
plain corresponds to <ft = 0. The hard partonic cross section (defined here for the gluon-gluon 
scattering) is of the form 
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£ cos(m0)/ m (AF,p) (2) 
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^ = log(^), P = log[3), (3) 



where a s = 3a s /n, 

AY = W ( 

k\k 2 

and I m is the solution of the BFKL equation with the consistency condition imposed. In this 
paper we shall use a fixed strong coupling constant corresponding to a s = 0.15. The solution 
may be written 

I m (AY,p) = / du / %L -e" AY cos(pv) (4) 
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where the kernels Xm are given by[] 

Xm(w,i/) = fi. 



(5) 



2 y r V 2 

with if} as usual, the logarithmic derivative of the Euler Gamma function. 

We evaluate the dw integral in (flj) by picking up the leading singularity uj^iy) which is the 
solution of the implicit equation 



(cf.(^D), to obtain 



I m (AY, p) = / dvRl(v) exp(u» m (v) AF) cos(H 
Jo 



where the residue R is 
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dXm(uJ,ls) 
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Hence at LO the above solution @ would simplify to 

wlM = Xm(0,f), 



with = 1. 
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Figure 1: uj^u), for < m < 4. In decreasing order, the curves correspond to m = 0, 1, ... 4. 



In Fig. [1] we illustrate the solution 00^(1/), (§), for various values of m. For AF ^> 1 the 
dominant contribution comes from the region v ~ 0. If we then expand around the saddle 
point, v = 0, we see that 

UAy,p) ~ exp(^(0)AF). (10) 

x An alternative resummation of subleading log(l/a;) effects, based on keeping the collinearly enhanced con- 
tributions [ fl3| is found to give the same result if ip(z) of (|^) is approximated by 1/z, apart from the inclusion 
in pa] of running a s and non-singular terms in the splitting function. 
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The total forward+backward jet cross section is controlled by the m = contribution and 
the intercept Wq(0) gives the rate of growth as the rapidity interval increases. We emphasize 
that the "subleading" intercept uJq(0) is significantly smaller than its LO counterpart, and 
significantly above the NLO approximation. 

It can also be seen from Fig. |l| that the intercepts ^(0) decrease with increasing m. This 
implies the dominance of the m = term at large values of AY, which determines the rate of 
decorrelation of the azimuthal angle between the forward and backward jets. In practice we do 
not restrict ourselves to the saddle point approximation, but rather evaluate (0) numerically. 
In Fig. 2 we plot the azimuthal distributions da/d<f) 

-£ = J dx 1 dx 2 dkldkl ——^— 2 (11) 

where the integration region A is restricted by the experimental cuts j| . To be precise the cuts 
which we impose are as follows: 

(i) the tagged jet tranverse momenta: |fej| > 20 GeV, and max(|fci|, |fc 2 |) > 50 GeV; 

(ii) the tagged jet pseudorapidities \rji\ < 3.0 and the pseudorapidity distance between jets 
|Ary = \i]2 — was required to fall into one of two windows: 2.5 < At] < 3.5 or 
4.5 < Ar] < 5.5. 



Fig. 2 shows the results corresponding both to the LO BFKL equation and those which are 
obtained from the BFKL equation with the subleading effects generated by the consistency 
constraint. We confront our predictions with the experimental data obtained at the Tevatron. 
We see from these figures that the subleading effects reduce azimuthal decorrelation obtained 
from the LO BFKL equation. The resulting distributions are however still too flat when 
compared with the experimental data. 

It has been pointed out in ref. Jl5| that the approximations which are adopted within the 



BFKL formalism may not be accurate, particularly in the region away from the back-to-back 
configuration of the two jets. It has been proposed that in order to improve accuracy of the 
predictions one should use the effective rapidity difference y(m) for each azimuthal projection 
m in the place of Ay in the solutions I m of the BFKL equation(s). The effective rapidity 
intervals y(m) are defined as: 

y(m) - AY - — cos ( m< ^ ( dc r l dyidy2dk 1 dk 2 d(f>) 
J d(j> cos(m(j))(dao/ dyidy2dkidk 2 d(f)) 

where da / dy±dy2dkidk2d(f) and do§j dy\dy2dk\dk2d$ denote the exact 0(a^) contribution to the 
three jet cross section and its BFKL (i.e. large Ay) approximation respectively. 

Guided by this proposal we have used as the asymptotic variable the effective rapidity 

mr, *?, 4 0) = ^ T'lTT^fTl < 13 > 

da /{dyidy2dkidk2d(p) 
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Figure 2: Azimuthal decorrelation of the forward and backward jets for two rapidity separation 
intervals: (a) 2.5 < Ai] < 3.5, (b) 4.5 < At] < 5.5. The DO data Q are shown by points with 
errorbars. The BFKL (LO) and BFKL (CC) predictions are shown by triangles and circles 
respectively. 
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where da / (dyidy 2 dk 1 dk 2 d(p) should now contain the subleading BFKL effects. They just cor 
respond to the replacement: 

1 



AY -> Ay 

2 




(14) 



in the LO BFKL expression for do§j (dy\dy2dk\dk 2 d<$) . 

The use of effective rapidity y(AY, kf, <fi) defined by equation (p~3| ) aims at correcting 
the large AY approximation of the 2 — > 3 jets cross-section integrated over the rapidity of the 
third jet. In this approximation the integrated 2^3 jets cross-section is just proportional 
to AF. This is the result of the fact that the rapidity y^ of the third jet varies between the 
rapidities y\ and y 2 of the tagged jets (AY = 2/2 — and that, in the large Ay approximation, 
potential dependence of the momenta of the colliding partons on y 3 is neglected. The analysis 



performed in ref. [JOd] shows that this is an inaccurate approximation and that the use of 
the exact kinematics introduces a large suppression. The effective rapidity y(AY, k\ , k%, 4>) 
incorporates this suppression and reflects the effective rapidity range spanned by the (mini)jets 
radiated with rapidities between y% and yi In the formulae defining the exact 3-jet cross 
section a fixed coupling was used corresponding to the value of a s which entered the BFKL 
kernels. Let us mention, that there is a significant dependence of y(AY, kf , k%, <p) on the (f) 
variable. The conditional maximum with other variables fixed is always reached for </> = and 
if k\ ~ k\ then y(AY^k\^k\^ = 0) ~ Ay. We have noticed that the numerical convergence 
of the integral in (0) is poor if y is below 1.2. In this region we have therefore used the 
perturbative expansion of the solution of the (modified) BFKL equation up to 0(a z s ). This 
required a careful treatment of the 3 jet cross section for the third jet with low transverse 
momentum fc 3 . Integration over the third jet phase space leads to infrared divergence which is 
cancelled by the virtual correction to the Born term, when integrated over the region of phase 
space containing the singular point k% = 0. 

The results of the calculations incorporating the effective rapidity are summarised in Fig. 3. 
These plots show the distributions which are obtained from the BFKL framework with the 
subleading effects taken into account and compare with effective rapidity y(AY, kf, fcf, <fi) in 
the place of Ay as an asymptotic variable. Again we compare with the available experimental 
data. We see that the use of the effective rapidity significantly improves agreement with the 
data. In Figures 3(a,b) we also show results obtained from the LO BFKL equation using the 
corresponding effective rapidity. In contrast to the results of the BFKL calculations presented 
in Figures 2 (a,b), when the effective rapidity is used the azimuthal decorrelation is far from 
being saturated (which would correspond to the distribution flat in (f>) and it should be more 
sensitive to the particular form of the BFKL kernel. Thus, it is surprising that these LO BFKL 
predictions do not differ appreciably from those obtained using the BFKL formalism with 
subleading effects and the effective rapidity. This implies that the implementation of effective 
rapidity suppresses the potentially large effect of subleading corrections. It is mainly caused by 
the fact that the absolute values of the effective rapidities at the tails of the distributions are 
rather small (i.e. y ~ 1), where the cross-section is dominated by the 2 jet + 3 jet production. In 
fact, the difference which is visible at the tails of the distributions obtained from the BFKL LO 
and BFKL CC equations does not come from the difference in the cross sections at the tails 
but rather from the region of central which influences the tails through entering the common 
normalisation factor N. Finally, we have checked, that the sensitivity of the results to variations 
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Figure 3: As for Fig. 2 but with the effective rapidity y of ([13]) incorporated in the BFKL 
predictions. 
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Figure 4: The average value of cos0 for various At] intervals. The DO data (points with 
errorbars) compared with the BFKL LO calculations (triangles) and with the BFKL CC calcu- 
lations (circles). The band between the At] axis and the dotted line represents the correlated 
systematic error of the data. 

of the fixed a s between 0.13 and 0.18 is small and a s = 0.15 gives the best description of the 
data. 

With experimental cuts used in the DO analysis the cross-section for cental bins of </> contains 
contributions from configurations with ki + k 2 = 0. In order to describe properly these config- 
urations the full BFKL resummation is necessary. This is due to the fact that the kinematical 
suppression of additional gluon emissions is weak in this case, and the corresponding effective 
rapidity y is close to AY. Thus, a s y becomes large and it is not enough to retain only the first 
terms of the perturbative expansion. In other words, the 0(a>l) cross-section, calculated as the 
sum of the LO 3 jet and NLO 2 jet cross-section, is expected to be subject to a series of large 
(a s y) n corrections. In fact, we have checked explicitly that without the resummation, the data 
in the central <p bins are not well reproduced. 

We also studied the distribution of (cos0) in the bins of pseudorapidity separation of the 
jets A-q. Here we used only the method and parameters which gave the best description of the 
azimuthal angle distributions of the Mueller-Navelet jets. The results of our calculation are 
compared to the experimental data [|J in Fig. 4. Both the LO BFKL and CC BFKL equations 
with the effective rapidity y give a resonable description of the experimental data, when the 
statistical and the correlated systematical errors of the data are taken into account. However, 
it is clear that some small discrepancy between the theoretical and experimental values still 
remains. 

To summarise we emphasize the following results of our analysis: 
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(i) The BFKL formalism in its original form fails to describe the data on the azimuthal 
decorrelation of dijets. The dominant subleading effects do improve the results obtained 
within the leading logarithmic approximation, but they still generate more decorrelation 
in the tails of the azimuthal distribution than that required by the data. 

(ii) Most of the relevant contribution to the decorrelation in the Tevatron domain is present 
already in the fixed order 0(af) expressions. 

(iii) In order to describe the data it is essential to use the effective rapidity y defined by 
eq. ( |i"3D and to use the 2 jet (LO + NLO) + 3 jet cross-section in the domain of low 
values of y. In this region the BFKL formalism does not therefore introduce any new 
information, which would not be present in the fixed order calculation. 

(iv) The only places where the use of the complete BFKL framework seems to be necessary 
are in the very central bins of 0. 

(v) We therefore conclude that it is possible to describe the decorrelation data using the fixed 
order perturbative expressions in the tails of the azimuthal distribution and the complete 
BFKL resummation for <fi « 0. 
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